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ABSTRACT 

The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope discovered a rapid (~ 
5 days duration), high-energy (E > 100 MeV) gamma-ray outburst from a source identified with the blazar 
PKS 1502+106 (OR 103, S3 1502+10, z=1.839) starting on August 05, 2008 (- 23 UTC, MJD 54683.95), 
and followed by bright and variable flux over the next few months. Results on the gamma-ray localization and 
identification, as well as spectral and temporal behavior during the first months of the Fermi all-sky survey are 
reported here in conjunction with a multi-waveband characterization as a result of one of the first Fermi multi- 
frequency campaigns. The campaign included a Swift ToO (followed up by 16-day observations on August 07- 
22, MJD 54685-54700), VLBA (within the MOJAVE program), Owens Vafley (OVRO) 40m, Effelsberg-lOOm, 
Metsahovi-14m, RATAN-600 and Kanata-Hiroshima radio/optical observations. Results from the analysis of 
archival observations by INTEGRAL, XMM-Newton and Spitzer space telescopes are reported for a more 
complete picture of this new gamma-ray blazar. PKS 1502+106 is a sub-GeV peaked, powerful flat spectrum 
radio quasar (luminosity at > 100 MeV, L^, is about 1.1 x 10*^ erg s^^, and black hole mass likely close to 
10^ Mq), exhibiting marked gamma-ray bolometric dominance, in particular during the asymmetric outburst 
{L^/Lopt ~ 100, and 5-day averaged flux ^EywoMcW — 2.91 ± 1.4 x 10^^ ph cm~ s~^), which was 
characterized by a factor greater than 3 of flux increase in less than 12 hours. The outburst was observed 
simultaneously from optical to X-ray bands (Fo.3-10 kcV — 2.18lg ^ 10^^^ erg cm^^ s^^, and hard photon 
index ^ 1.5, similar to past values) with a flux increase of less than one order of magnitude with respect to 
past observations, and was likely controlled by Comptonization of external-jet photons produced in the broad 
line region (BLR) in the gamma-ray band. No evidence of a possible blue bump signature was observed 
in the optical-UV continuum spectrum, while some hints for a possible 4-day time-lag with respect to the 
gamma-ray flare were found. Nonetheless, the properties of PKS 1502+106 and the strict optical/UV, X- and 
gamma-ray cross correlations suggest the contribution of the synchrotron self Compton (SSC), in-jet, process 
should dominate from radio to X-rays. This mechanism may also be responsible for the consistent gamma-ray 
variability observed by the LAT on longer timescales, after the ignition of activity at these energies provided 
by the BLR-dissipated outburst. Modulations and subsequent minor, rapid flare events were detected, with a 
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general fluctuation mode between pink-noise and a random-walk. The averaged gamma-ray spectrum showed 
a deviation from a simple power-law, and can be described by a log-parabola curved model peaking around 0.4- 
0.5 GeV. The maximum energy of photons detected from the source in the first four months of LAT observations 
was 15.8 GeV, with no significant consequences on extragalactic background light predictions. A possible radio 
counterpart of the gamma-ray outburst can be assumed only if a delay of more than 3 months is considered 
on the basis of opacity effects at cm and longer wavelengths. The rotation of the electric vector position angle 
observed by VLBA from 2007 to 2008 could represent a slow field ordering and alignment with respect to the 
jet axis, likely a precursor feature of the ejection of a superluminal radio knot and the high-energy outburst. 
This observing campaign provides more insight into the connection between MeV-GeV fliires and the moving, 
polarized structures observed by the VLBI. 

Subject headings: gamma-iays: observations - quasars: individual: PKS 1502H-106- quasars: general - galax- 
ies: active - galaxies: jets - X-rays: galaxies 
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PKS 1502+106: a new and distant gamma-ray blazar in outburst discovered by the Fermi LAT 
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1. INTRODUCTION 

The Large Area Telescope (LAT), on boar d the Ferm i 
Gamma-ray Space Telescope (formerly GLAST; Ritz'' 2007l) . 
was successfully launched by NASA on 2008, June 11, from 
Cape Canaveral, Florida, on a Delta II Heavy launch vehi- 
cle. While still in the commissioning and checkout phase, 
it discovered and monitored bright, flaring gamma-ray emis- 
sion above 100 MeV from a source identified with the blazar 
PKS 1502+106 (historically also known as OR 103 and S3 
1502+10). The large field of view, effective area and sensitiv- 
ity and the nominal survey observational mode make Fermi- 
LAT an unprecedented all - sky monitor of 7-rav flares and 
source variability (see, e.g. 'McEne ryl2006tlThompson 120061: 
iLott et al. 2007; Atwood et al. 2005^^ 

At the beginning of August 2008, PKS 1502+106 was the 
second brightest extragalactic source in the 7-ray sky, ex- 
hibiting a sudden high-energy outburst announced in ATel 
#1650. This outburst successfully triggered the first (un- 
planned) Fermi multi-frequency campaign. Major renewed 
gamma-ray activity observed by Fermi in January 2009 was 
announced via ATel #1905. 

PKS 1502+106 is a luminous, quasar-like (optically broad- 
line and flat radio spectrum) AGN discovered during the 178 
MHz pencil beam survey from the Mullard Radio Astronomy 
Observatory, Cambridge, UK, (appearing in a list not included 
in the 4C catalog; Crowther & Clarke 1966; Williams et al. 
IT967h . and was re-observed and characterized as an extra- 
galactic source by both the Australian Nation al Radio Astron- 
omy Observatory of Parkes, NSW, Australia, dPav et alj 19661 
id.: PKS 1502+106), and the Ohio State Universit y ("Big 
Ear") Radio Observatory, Delaware, OH, USA, (Fitc h et al] 
[1969, id.: OR 103). The source exhibited substantial radio 
flux variations (factor > 2), a high degree of linear polariza- 
tion, a core-dominated, one-sided and curved radi o jet with 
a variable, a com plex morphology at VLBI scales (l An et alj 
120041: iLister et a l. 2009a), and a compact large scale struc- 
ture. 1 1 VLBA observations at 15.4 GHz performed between 
Aug. 1997 and Aug. 2007 showed a FWHM major beam 
axis in the range 1.02-1.57 mas, a minor axis beam axis of 
0.5 mas , and a total flux density in the range 0.88-1.93 Jy 
( Lister e t al. 2009a), and apparent jet speed of (14.8 ± 1.2)c 
( iLister et aLi2009 b). The 22 and 37 GHz flux history shows 
several long-term flares (> 1 Jy variations, i.e. ^ 60% of the 
total flux range span, on typical timescales of a year, and peak 
fluxes well above 2 Jy), with at least five flares and an aver- 
age trend that was sligh tly increasing from 1988 to mid-2004 
jTerasranta et alj|2005h . WMAP fluxes at similar frequencies 
(K, Ka, Q bands) are in agreement with these flux ranges 
jLopez-Caniego et al.l 120071) . The Doppler factor estimated 
from t he observed 37 GHz variability and brightness temper- 
ature (jHovatta e Tal. 2009) agrees with the jet speed (14.6c) 
cited above, a Doppler factor Vyar = 12 and viewing angle 
— 4 7° 
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This radio blazar was identified in th e optical band by 
iBlake (I1970R with a position refinement by I Argue & SuUivanl 
(119801) . while an initial spectrosc opic in spection was per- 
formed by [Burbidge & Strittmatter (119721) . Variations > 2.5 
mag were observed in its optical flux history (Palomar-Quest 
and Catalina Sky Surveys, ATel #1661), together with a vari- 
able and relatively high degree (up to 20%) of linear polar- 
ization, pointing out a dominant synchrotron emission with 
no observed dilution by thermal components. The redshift of 
PKS 1 502+ 1 06, as confirmed by the good S/N spectrum of the 
Sloan Digital Sky Survey (z = 1.8385 ± 0.0024 at high con- 
fidence), is in agreement with the value z = 1.839 estimated 
previously bv lSmith et alj (Il977h . A less remote value (z = 
0.56) is reported in oth e r works dBurbidge & Strittmatteii 
[T97a iWright etal.l[T97l iWilkes et al.lll983 t). although the 
possible multiple Mgll absorption system (pointed out by a 
feature shortward of the 4388 A emission line) would be very 
unusual for a low redshift object. 

Serendipitous X-ray data of PKS 1502+106 are available 
because the source lies about 7' NE of the bright Seyfert type- 
1 galaxy Mkn 841, althou gh only one multifrequency work 
dedicated to this blazar (George et al. 1994) has appeared pre- 
viously. Early X-ray observations (ROSAT, ASCA) showed 
low-amplitude variations on short timescales (factor > 2 on 
timescales of a year), a flat 0.1-10 keV photon index Tx 
between 1.4 and 1.9, and an intrinsic X-ray luminosity of 
i2-iokov = 1-2 X 10''^ erg s"\ and a 2-10 keV flux in the 
range 4.9 - 6.54 x 10~^^ erg cm~^ s~^G eorge et al. 1994 
lAkivama et afll200l Imtanabe et alJl2004l) . PKS 1502+106 
was speculated to be a possible 7-ray source before the 
LAT discovery because of the super luminal motions of jet 
components (u p to 187 ± 15/ias/year lLister & Homanll2005t 
lAn et al.ll2004l) . and the multiwaveband spectral indexes arx 
and apx (co nsistent with other FSRQs detected by EGRET, 
iGeorgeetai. 1994). Only modest intrinsic X-ray absorption 
was suggested by this work, and the optic al and near-IR red- 
dening claimed in IWatanabe et al.l (120041) is probably due to 
the synchrotron jet dominance at these low frequencies rather 
than by absorption from inner nuclear light. 

A relation involving the misalignment between the pc- 
and kpc-scale radio structure (positio n angle) and the 7 - 
ray emission was postulated as well (Cooper et al.l 120071) . 
However, only a cumulative 2a upper limit by EGRET 
of 7 X 10~* ph cm~^ s^^ was reported (Phase/Cycle I, 
combi ned Viewing Perio ds: 24.0 to 25.0, i.e. April 02-23, 
1992; Fichtel etalj |1994), and the source was likewise unde- 
tected in the following EGR ET cycles ( Hartman et alJ[T999l ; 
ICasandjian & Grenieill2008h . 

In the following we use a ACDM (concordance) cos- 
mology with values g iven within Icr of the WMAP results 
dKomatsu et al.l 120091) . namely h = 0.71, = 0.27 and 
J7a = 0.73, and a Hubble constant value Hq = lOOh km s^^ 
Mpc-i. 

In Section 2, first results on the 7-ray identification, the ob- 
served MeV-GeV outburst and the subsequent four months 
of monitoring by the ferm/-LAT are described. In Sec- 
tion 3, multifrequency results obtained through simultaneous 
optical-UV-X-ray observations by Swift (thanks to a 16-day 
long monitoring following a triggered Target of Opportunity, 
ToO), and by radio-optical observatories (the 40m dish tele- 
scope of the Owens Valley Radio Observatory, the Effelsberg 
100m dish radio telescope, the ring radio telescope RATAN- 
600, the VLBA within the MOJAVE program, and the Kanata 
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Fig. 1. — Left panel: LAT count map cumulated on a nine-month (Aug. 
2008 - Apr. 2009) baseline, weighted and smoothed by the point spread 
function (PSF) such that higher energy photons are mapped to higher in- 
tensities. The map is in arbitrary units in the energy range 0.1-100 GeV and 
in a 2.5° X 2.5° region centered on PKS 1502+106. The qualitative circle 
sizes of the PSF at 200MeV, 2GeV and 20GeV are outlined for reference. 
Right panel: LAT source localization with 95% and 68% uncertainty radii 
(red circles) superimposed on an arcmin-scale (R-band) optical image show- 
ing also the X-ray counterpart eiTor box by the Swift-XKI observations and 
the radio position and intensity contours by VLA of PKS 1502-1-106. The 
best LAT source position, calculated on the same nine-month period, with 
the pointlike tool is RA: 226.10179°, Dec: -1-10.4927°, A = 0.0027°, 
with 68% and 95% LAT error circles of 0.0077° and 0.0124° respectively. 

telescope of the Higashi-Hiroshima Observatory) are summa- 
rized. In addition, past and unpublished observations by the 
XMM-Newton and Spitzer space telescopes are analyzed and 
presented in Section 4 for a more complete picture. Finally, 
in Section 5 and 6, discussion and concluding remarks are re- 
ported. 

2. GAMMA-RAY OBSERVATIONS AND RESULTS BY FERMI-LAT 
2.1. LAT observations 

The LAT instrument is a pair tracker-converter telescope 
comprising a modular array of 16 towers — each with a 
tracker based on silicon micro-strip detector technology — 
and a calorimeter based on a hodoscopic array of 96 CsI(Tl) 
crystals, surrounded by an Anti-Coincidence Detector capa- 
ble of measuring the directions and energies of cosmic 7- 
ray photons with energies from 20 MeV to > 300 GeV (for 
details, see, e.g. [Bellazzini et al. 2002; Michelson 200'^ 
lAtwood et al.ll2007ll2009l:IAbdo et al. 2009h). 

The reduction and analysis of LAT data was performed us- 
ing the Science Tools v.9.8, based in particular on a unbinned 
maximum-likelihood estimator of the spectral model param- 
eters (gtlike tool). Events were selected using the Instru- 
ment Response Functions (IRFs) P6_V1 .DIFFUSE. This se- 
lection provides the cleanest set of events (in terms of direc- 
tion, energy reconstruction and background rejection) at the 
cost of reduced effective area at low energies, and takes into 
account the differences between front- and back-converting 
events. To minimize contamination by Earth albedo 7-ray 
events that have reconstructed directions with angles with re- 
spect to the local zenith > 105° have been excluded. For this 
object with high Galactic latitude, events are extracted within 
a 10° acceptance cone centered at the PKS 1502H-106 radio 
position. This cone, substantially larger than the 68% contain- 
ment angle of the PSF at the lowest energies, provides suffi- 
cient events to accurately constrain the diffuse emission com- 
ponents. The gtlike model includes the PKS 1502H-106 
point source component, two other point sources from the 3 
month catalog (both faint and low-confidence sources with 
TS ~ 0.9% of the TS value of PKS 1502H-106 for the same 
period), a component for the Galactic diffuse emission (GAL- 
PROP code, see, e.g. .Moskalenko et al..2003. and references 



therein), and an isotropic component including the extragalac- 
tic diffuse emission and the residual background from cosmic 
rays. 

The Fermi-LAil data of PKS 1 502h- 1 06 presented here were 
obtained during the first four months of the LAT survey (Aug.- 
Dec. 2008). In this period PKS 1502H-106 was one of the 
most persistently bright, variable sources in the high-energy 
sky and almost certainly the source with the highest luminos- 
ity. The background contribution within a few degrees was 
only a small fraction of the source count rate, with no nearby 
source confusion. The time interval was sufficient for a fine 
determination of the average spectrum, for a first look at the 
mid-timescale variability and detection of posterior flares, for 
a refined localization, and a first cross comparison with the 
other multifrequency monitoring data. The first 7-ray detec- 
tion of PKS 1502H-106 by the LAT occuiTed in July 2008, 
when it was confirmed by the high-level Automatic Science 
Proce ssing pipeline monitoring (ASP; iChiang et al.l 120061 
l2007h . based on a wavelet-based (pgwave) quick-look de- 
tection tool (e.g. iDamiani et al.llT997t iMarcucci et al.ll2004t 
ICiprini et al.ll2007al) and a maximum likelihood analysis, and 
by the LAT Source Catalog algorithm (mr_filter) based on 
wavelet analysis in the Poisson regime (Starck & Pierre 199^ 
and t he peak-finding tool sExtractor (iBertin & ArnoutsI 
119961) . The rapid and markedly time-asymmetric 7-ray out- 
burst announced in ATel #1650 and triggering an unplanned 
Target of Opportunity (ToO) multifrequency campaign was 
seen from Aug. 05 (about 23 UTC) until about Aug. 11, 2008 
(^ 5 days duration), with a fast rise, slower decay, and an 
approximately two-day sustained peak flux. 

Some caveats related to the preflight instrument response 
functions (P6_V1), which overes timated the acceptan ce at low 
energies, are briefly described in I Abdo et al.l (l2009bl) . 

2.2. Gamma-ray source localization, association and 
identification 

The LAT PSF and sensitivity provides an unprecedented 
angular resolution in gamma-rays (68% containment ra- 
dius better flian ^ 1° at 1 GeV, [Atwood et al. 200"^ 
ICecchi et al.ll2007l; lAtwood et al.ll2009l; lAbdo et al. 2009c h), 
making the association and identification processes less dif- 
ficult than in previous experiments. In the case of this very 
bright 7-ray source, we obtained — beyond the good spatial 
"association" — a firm "identification" with PKS 1502 -H106. 
The 3-month bright source list results (" Abdo et ani2009 c, id: 
OFGL J1504.4H-1030), provided a good initial localization: 

RA: 226.12°, Dec: H-10.51°; rgs = 0.05° and VtS = 88.2, 
(rgs being the radius of 95% confidence region, TS the like- 
lihood test statistic from the 200 MeV to 1 00 GeV analy- 
sis). A pplication of the point like tool (iBurnettI 120071: 
lAbdo e t al. 2009c) on a much longer (nine-month, Aug. 2008 
- Apr. 2009) LAT dataset with very high statistics, provided 
an excellent localization (outlined in Fig[T|: RA: 226. 10179°, 
Dec: H-10.4927°, A 0.0027°, with 68% and 95% LAT error 
circles of 0.0077° and 0.0124° respectively (statistical only). 
Studies of bright source localizations indicate a systematic un- 
certainty in the localization of < 30", that can be taken as an 
estimate of the systematics with this tool (more details on the 
production o f pointlike den sity maps and localization are 
described in ICamilo et al. 112009) . The relevant improvement 
in the localization carried out on a nine-month baseline is due 
to the high variability that occurred with this source. These 
localization values are in agreement with the VLBI radio and 
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Fig. 2. — Main panel: Likelihood flux (E>100MeV) light curve obtained in daily bins from Aug. 02 to Dec. 15, 2008. The outburst state and the subsequent 
post-flare (a lower and intermediate level brightness, far from the faintest state observed) period with simultaneous monitoring by Swift are represented by the 
two horizontal lines. Left inset plot: A zoom on the corresponding flux light curve around the outburst period obtained using finer, 12-hour bins (lower statistics). 
Right inset plot: The gamma-ray (E>100MeV) photon index values for the same period using daily bins as in the main panel light curve. 



optical positions of PKS 1502-1-106, the VLA contours and 
the Swift XRT error box (Fig. [IJ. PKS 1502H-106 is the only 
bright VLA radio source (calibrator source list) located within 
the LAT 95% confidence circle. The Seyfert galaxy Mkn 841 
(observed with a hard X-ray cutoff, see SectlDl is positioned 
well outside of these localization circles, as are other HB89- 
catalog quasars in this region. 

Beyond the excellent spatial association, the most secure 
and distinctive signature for firm identification of this new 
gamma-ray source found by Fermi is the observed correla- 
tion between the 7-ray, X-ray and optical-UV variability (see 
sections |3] and |4]i. This object was also a member of the 
pre-launch CGRaBS (iHealey e t al. 2008 . object id: CG RaBS 
J1504H-1029) and Roma-BZCAT (M assaro et al.ll2009l object 
id: BZQ J 1504-1- 1029 ) catalogs listing candidate gamma-ray 
blazars. Fi nally, a method based on a "figure o f merit" (de- 
scribed in ISowards-Emmerd et al.ir2003l l2005h for this LAT 
source position provides a very high Ukelihood of identifica- 
tion with PKS 1502H-106. 

2.3. Gamma-ray temporal behavior 

The typically bright 7-ray flux and the enduring activity 
shown by PKS 1502-1-106 in 7-rays, allowed a firm detec- 
tion of the source by the LAT on a daily basis. Fig. |2] 
shows the light curve (daily bins, E > lOOMeV) extracted 
with the gtlike tool over the first four and half months 
of LAT all-sky survey. A fast-rising, markedly asymmet- 
ric and bright outburst was found, with a factor > 3 of in- 
crease in flux in less than 12 hours. The integrated flux at 
E > 100 MeV averaged in the 12h bin of the peak emission 
(between Aug. 05 and Aug. 06, 2008, i.e. DOY 218-219) 
was Fe>ioo MoV = (3.7 ± 0.7) x lO"*^ ph cm"^ g-i (statis- 
tical only), as measured when the LAT instrument was still in 
commissioning and checkout phase (all-sky nominal mode). 
The emission from PKS 1502-1-106 then faded more slowly in 
the following days, and the entire outburst interval spanned 
Aug. 05 around 23 UTC until to Aug. 1 1 around 00 UTC, 
2008 (i.e. DOY, 218.95-224.0, - 5 days duration Fig. EJ. 
The peak flux appeared elevated for less than two days, rival- 
ing the brighest a ppar ent flux of other extragalactic sources at 
that time (Section uTt . The finer, 12-hourbin light curve (^ 8 
Fermi orbits, ensuring still similar exposures) reported in the 
right inset panel of Fig. |2] shows the elevated flux held for 
at least one 1 .5 days, while the slower fading decay exhibits 
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Fig. 3. — Time series analysis of the LAT light curve presented in Fig|2] 
the periodogram, first order structure function (upper inset) and power spec- 
trum (lower inset). These functions show a variability with a power spectrum 
consistent with 1//^ '* fluctuations. This indicates a variability mode placed 
between flickering and shot noise. The horizontal dashed line represents the 
0.01 false alarm probability threshold (99% significance that the signal de- 
tection is not wrong). 

a high flux (-Fe>ioo MeV > 2 x 10~^ ph cm~^ s^'^) plateau, 
or a secondary superposed flare, that extended for about 2.5 
days. During this outburst a "harder when brighter" spectral 
trend is suggested, despite the photon index error dispersion 
(see the daily photon indices reported on the right inset panel 
of Fig. |2]l. The gtlike performances and the current IRFs 
used did not allow to go below a daily binning in the extrap- 
olation of the photon index values, and this was possible in 
an acceptable way with respect to amplitude of the statistical 
error only for the high-flux and high photon count statistics 
available around the outburst epoch. 

A consistent level of variability, with a couple of minor but 
relevant rapid flares, occurred after the major outburst with 
fluctuations on timescales of weeks. Renewed activity and in- 
creased average brightness from the end of Nov. 2008 were 
observed. Two rapid flares approached a maximum peak flux 
of i^£;>ioo MoV ^ 2 xlO^^ ph cm^^ s^^(daily bin estima- 
tions), on Sept. 06 (DOY 250), where a simultaneous cross- 
correlated optical flare was observed as well, and on Oct. 30 
(DOY 304). Visual inspection of the light curve reported in 
Fig|2] suggests a period of higher activity beginning in mid 
Nov. 2008 (after about DOY 320), and, in general, a series of 
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modulations occurring on about a one-month timescale, with 
faster fluctuations and rapid flare episodes superposed. The 
power spectral density (PSD) shows a power-law dependence, 
P{f) oc Similar time scale dependence is exhibited 

by the first order structure function (Fig. [31 and by the au- 
tocorrelation function. More detailed variability analysis for 
th is and other blazars using a longer dataset will be presented 
in lAbdo et al](l2009dh . 

2.4. Gamma-ray spectra 

We have analyzed the time-averaged spectra of PKS 
1502H-106 for three intervals: the high-state of the outburst 
(DOY, 218.95-224.0, i.e. Aug. 05 - Aug. 11, 2008, about 
5 days); the post-flare period characterized by an intermedi- 
ate brightness level and during which simultaneous Swift ob- 
servations were performed (DOY 224.0-235.42, i.e. Aug. 11- 
Aug. 22, 2008, about 11.4 days); the longer and heteroge- 
neous period that includes the outburst and the following 4 
months of Ferm/-LAT monitoring (Aug. 02 to Dec. 07, 2008, 
^ 126 days, where the source displayed different stages of 
activity and significant variability). Events below 200 MeV 
were excluded from these analyses because of calibration un- 
certainties at those energies. An isotropic background model 
was used as PKS 1502H-106 was very bright relative to other 
point sources during the period stated above, because is lo- 
cated at high Galactic latitude, and because checks with more 
complex models provided no significant difference. Further- 
more, no appreciable differences were observed using differ- 
ent acceptance cone radii for the event extraction. 

The spectra for the post-flare and cumulative 4-month 
datasets can be consistent with a log-parabola (LP) model, 
dN/dE (X ^-("+/3i°g ( ^)) (see, e.g . Landau etal. 198 
Inoue & Tak aharal 11996^ |Fossati etall 12000; Massaro et al. 



2004t iPerlm an et al. 20^^ The likelihood ratio test 
jMattox etal.iil996i) rejects the hypothesis that the spectrum 
is a power law (null hypothesis) against the one that the spec- 
trum is curved as a log-parabola model. This model is pre- 
ferred over a simple power-law model at the 1 1 .4 sigma sig- 
nificance level. Broken power-law (BPL) fits show a similar 
improvement over single power-law models, but we find no 
evidence to prefer BPL over the log-parabola representation. 
For the full time interval characterized by very high statistics, 
the logarithm of the likelihood increases significantly when 
allowing (3 to vary, and an increase of the value for the BPL 
with respect to simple power-law model of the same order of 
the increase for the LP vs. power-law test is not observed. 



TABLE 1 

Summary of the unbinned likelihood spectral fit 
ABOVE 100 MeV 



Interval [MJD (DOY)] Best-fit Model and Parameters 



All observations 
54682.680 (217.680) 
54775.580(310.580) 



Log-parabola 
a = 1.94 ±0.05 
13 = 0.10 ±0.02 

FiS>ioo McV = 6.90 ± 0.34 x 10"'^ [ph cm"^ s"^ 



Outburst/high state 
54683.955 (218.955) 
54688.985 (223.985) 



Power Law 

r = -2.06 ±0.017 

Fis>ioo McV = 29.1 ± 1.4 X 10-^ [ph cm-^ i 



Post-flare/lower state 
54689.063 (224.063) 
54775.580(310.580) 



Log-parabola 
a = 1.87 ±0.20 
/3 = 0.18 ±0.08 

Fis>ioo McV = 5.32 ± 1.03 x 10~'^ [ph cm'^ ; 



The LP description introduces the advantage of only one ex- 
tra parameter (while BPL model adds two parameters) with 
respect to the simple power-law model, it allows modeling 
of mild spectral curvatures with no abrupt cutoffs, and pro- 
vides a better phenomenological physical interpretation. On 
the other hand there can be still caveats when using gtlike 
with a broken power-law model, in particular in determining 
the break energy when statistics do not allow a high number 
of energy bins. It is also plausible that an energy spectrum 
averaged over a long period of time, and containing different 
activity stages with time varying hardness, may exhibit an ap- 
parent curvature. Finally this does not exclude BPL model if 
the spectrum is extracted in different time intervals. 

The average spectrum during the outburst state is consistent 
with a simple power-law model, dN/dE cx E^ . The outburst 
state shows a rather hard spectrum, suggesting a maximum 
peak in the MeV energy bands (in the vF^, representation), 
in agreement with the LP peaks found for the spectra cited 
above. The extrapolated and averaged fluxes integrated above 
100 MeV, and the spectral fit parameters for all three periods 
are shown in Table [T] 

3. SIMULTANEOUS MULTIFREQUENCY OBSERVATIONS 

Because of the reasonably uniform exposure and high sen- 
sitivity of the LAT, and the broad-band radio-to-gamma-ray 
emission of this kind of AGN, simultaneous multifrequency 
data are very important to the investigation of the physical 
properties of supermassive black holes and relativistic jets, 
beyond the benefit of a firm source identification (Section 
12.21 ). With this in mind, several campaigns on a few selected 
objects, or ToO list of candidates for flaring so urces, were 
prepared pre-launch by the Fermi collaboration (iTostil 120071 : 
iThompson 2007). PKS 1502H-106 was a previously unknown 
7-ray source, with no pre-planned multifrequency campaign. 
But following the LAT outburst (reported in ATel #1650), a 
ToO campaign was initiated on Aug. 07, 2008. This was 
the first Fermi multifrequency campaign that had not been 
planned pre-launch, and saw triggers for ToO pointings by IN- 
TEGRAL and Swift, long-term radio flux and structure moni- 
toring, as well as optical observations by ground based facili- 
ties. 

The fast response ToO pointing by INTEGRAL provided 
200ks of observations during the period Aug. 09, 01:53 UT - 
Aug. 11, 15:12 UT, 2008 (revolution 711). However, PKS 
1502H-106 was not detected (preliminarily) by the imager 
IBIS on board INTEGRAL. Extrapolating the X-ray flux ob- 
served by Swift, the hard-X-ray flux had likely already faded 
to slightly below the IBIS detection thresh old in this epoch 
(more details will appear in lPian et al.ll2009l) . 

3.1. Simultaneous X-ray and UV-optical observations and 
= results by Swift 

= The Swift satellite (iGehrels et all 12004 *) performed a ToO 
monitoring campaign of PKS 1502H-106 with daily snapshots 
from Aug. 07 to Aug. 22, 2008. This quite long-term ob- 

] serving campaign by Swift allowed extended daily snapshots 

Tor about 16 days, using the three instruments onboard: the 
X-ray telescope (XRT) for the 0.2-10 keV energy band, the 

_Ultraviolet/Optical Telescope (UVOT) for multiband photom- 
etry, and the Burst Alert Telescope (BAT) for the 15-150 keV 
hard X-ray band. BAT data were not used because of source 
confusion problems with Mkn 841 which is about a factor of 

LjO brighter than PKS 1502H-106 in the hard X-ray band. The 
16 days of observations by Swift allow for cross-correlation 
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Fig. 4. — Simultaneous gamma-ray and multifrequency light curves ob- 
tained during the multiwavelength campaign of August 2008 triggered by the 
high energy outburst discovered by Fermi-hAI. The flux above 100 MeV, 
the X-ray flux (0.3-lOkeV) by Swift-'KKT, the six-band fluxes monitored 
by Swift-WOT, the Kanata-TRISPEC differential relative magnitude hght 
curves (optical AV and near-IR A J bands) and con'esponding measures of 
the hnear polarization degree, and the 15 GHz radio light curve from OVRO 
40-m are reported. 

Studies between the 7-ray, X-ray and UV-optical bands during 
both the active flaring stage and the fading post-flare stage of 
PKS 1502H-106. 

The XRT was set in photon counting mode, and the data 
were processed by the xrtpipeline with the use of standard 
software (HE AD AS software, v6.4) and standard filtering and 
screening criteria. The XRT events in the 0.3-10 keV energy 
band were extracted using the XRTGRBLC FTOOL from cir- 
cular regions centered on the source position with variable 
radii depending on the source intensity and applying correc- 
tion for vignetting, Point-Spread Function corrections and hot 
pixels and columns with the use of exposure maps. The XRT 
X-ray flux light curve is shown in the second panel of Fig. 
H The 0.3-lOkeV count rate of PKS 1502H-106 measured by 
XRT was at a level 0.05 counts/sec (from our data), up from a 
level of 0.02 counts/sec (archival past observations). The 0.3- 
10 keV XRT 16-day long light curve obtained in August 2008 
(FiglU shows an initial count rate of 0.05 counts/sec, and a 
gradual decay down to the level of about 0.02 counts/sec. 
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Fig. 5.— Swift-XRT combined 0.3-10 KeV spectra of PKS 1502+106 ex- 
tracted for the high state (MJD: 54685-54689) and the subsequent low state 
(MJD 54690-54701), mapping the X-ray behavior simultaneous to the LAT 
flare and to the post-flare, relaxing activity and brightness. 



TABLE 2 

Analysis summary of the simultaneous data obtained by the 

XRT INSTRUMENT ON BOARD Swift. 



Obs. id. (date) 


Best-fit Model and Parameters 
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The X-ray spectrum of each observation segment was fit- 
ted with an absorbed power law. Because of the low num- 
ber of events from the source, event were not grouped and 
C-statistics was used, fixing the column density Nh to the 
Galactic value Nhi = 2.19 x 10^° cm~^ in that direction 
(in agreement with values u sed, for example, in iGeorge et alj 
Il994t lAkivama et alj|2003l) . and using z = 1.839. The error 
on the photon index and the flux (0.3-lOkeV) is large due to 
the low statistics. The background photons were selected in 
a circular region close to the source. No significant photon 
index variation was observed between the high and the low 
state, while the count rate and flux did vary. 

The Swift Ultraviolet/Optical Telescope (UVOT) photom- 
etry was done using the publicly available UVOT FTOOLS 
data reducti on suite, and is in the UVOT photometric system 
described in lPoole et al.l (|2008). The photometric data points 
were corrected for Galactic extinction using the dust maps of 
Schlegel et al. ( 1998) and the Milky Way extinction curve of 
iPeil (11992). These simultaneous multi-band optical and UV 
data show an increase of about 2 magnitudes in all filters when 
compared with the past-years archival values (i.e., from about 
19 to 17 in B band). The flux light curves in the six UVOT 
bands are shown in the third and fourth panel of Fig. |4] These 
fluxes appears to be well correlated. A slight rise in flux of 3 
days is observed in all the UVOT bands, followed by a fading 
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similar to the flux decrease seen in 7-rays and X-ray bands. If 
the time of the observed UV and optical maximum is related 
to the flare activity at higher energies, this would imply an 
interesting time lag of about 4 days. 

3.2. Simultaneous near-infrared and optical monitoring 

PKS 1502+106 was also monitored in the optical V and 
near-infrared J bands with some photometric and polarimetric 
snapshots by the TRISPEC istrument attached to the 1.5-m 
"KAN ATA" telescope at the Higashi-Hiroshima O bservatory, 
Japan ('Watana be et al.1 120051: lUemura et alJ 12008 1 within a 
twofold program of optical follow up for LAT flaring sources 
and regular monitoring of about 20 blazars. Imaging relative 
photometry was performed using some comparison stars in 
the same field, but due to the absence of an accurate calibra- 
tion for this field we prefer to report only the relative magni- 
tude difference Amag with respect to the minimum level (fifth 
panel Fig. |4|. The optical and near-infrared band imaging 
photometry is performed simultaneously in TRISPEC with a 
unit of polarimetric sequence (consisting of successive expo- 
sures at four position angles of the half-wave plate, where a 
set of linear polarization parameters, Q/I, U/I, are derived). 

These flux observations, performed for a longer interval 
with respect to the Fermi-Swift campaign (i.e. until Sept. 22, 
2008), show a high correlation between the V- and J-band 
light curves and show an optical decay phase comparable to 
that observed in the UVOT photometric observations (fifth 
panel of FiglH. Remarkably, a strong correlation with the 
LAT gamma-ray light curve is found, including the first (Sept. 
04-07, 2008) of the possible minor flares occurring after the 
initial large outburst. The observation of a flare in the opti- 
cal (V) and near-IR (J) bands, simultaneous with a second 7- 
ray flare have the twofold advantage of providing a validation 
of such a minor LAT flare as a real feature displayed by this 
blazar, and, even more crucially, in confirming the firm identi- 
fication of the new gamma-ray source seen by Fermi with the 
blazar PKS 1502+106. 

Comparing the Kanata-TRJSPEC V-band and J-band col- 
ors, the V — J color index varies between 2.05 and 1.69 
(during the Sept.04-07 minor, rapid flare cited above). On 
the other hand, the degree P of linear optical (in y, J bands) 
polarization observed (sixth panel of Fig|4]i , remains rather 
scattered by error dispersion irrespective the flux level, even 
during the minor flare mentioned (the maximum degrees 
recorded during the monitoring were P{V)max = 15 ± 3% 
and P(J),„ar. - 13 ±4%). 

3.3. Simultaneous radio flux-structure data by single-dish 
and VLBI observations 

As part of an ongoing blazar monitoring program, the 
Owens Valley Radio Observatory (OVRO) 40-m radio tele- 
scope has observed PKS 1502+106 at 15 GHz approximately 
every two days since mid-2007. Flux densities for the peri- 
ods from July 26 to September 3, 2008 (MJD 54673-5471 1) 
and October 23 to December 9, 2008 (MJD 54762-54809) 
are shown in Figure|6l Flux densities were measured using az- 
imuth double switching as described in Readhead et al.l (119891) 
after peaking up on-source. The relative flux density uncer- 
tainty for this source is dominated by a conservative 1.6% 
systematic error with a typical thermal error contribution of 
5 mJy. Abso l ute flu x density is calibrated to about 5% using 
the iBaars al.l (119771) model for 3C 286. This absolute uncer- 
tainty is not included in the plotted errors. 



The measured flux densities in the MJD 54673-54711 
(DoY 208-246) time period fit a 1.69 Jy constant-flux model 
with x^/{N - 1) = 0.70 {N = 15 data points). This in- 
dicates no statistically significant variability in this time pe- 
riod. The beginning of a bright radio flare is apparent in the 
MJD 54762-54809 (DoY 297-344) time period with an in- 
crease of at least 30% over the earlier mean flux density. 

A less intensive monitoring at 37 GHz was carried out with 
the 13.7m radio telescope at Metsahovi Radio Observatory, 
Helsinki University of Technology, Finland. The flux den- 
sity scale is set by observations of DR 21, and sources 3C 84 
and 3C 274 are used as secondary calibrators. A detailed de- 
scription on the data r eduction and analysis can be found in 
iTerSsranta et~ai] ( 119981) . 

The PKS 1502+106 flare was also followed-up by the Ef- 
felsberg 100-m radio telescope with four multi-frequency ra- 
dio spectra obtained on August 23, September 16, October 18, 
and December 6, 200 8 (within the F-GAMMA project, see 
iFuhrmann et al.ll2007t) . Each radio spectrum consists of si- 
multaneous measurements at various frequencies between 2.6 
and 42 GHz. The observations were performed using cross- 
scans in azimuth/elevation with the number of sub-scans 
matching the source's brightness at the given frequencies. The 
d ata reduction was done u sing standard procedure s described 
in IFuhrmann et all (l2008h : lAngelakis et all (120081) . 

Other radio observations are available at six frequencies be- 
tween 1 and 21.7 GHz thanks to the 600-m ring transit radio 
telescope of the Russian Ac ademy of Sciences RATAN-600 
(iKorolkov & Parijskijlll979h . which observed the source on 
September 10, 26, and October 2, 2008. A weighted average 
of these three observations, is presented in Fig |7] Previous 
RATAN-600 data which cover the period between 1997 and 
March 2008 are also shown for comparison. The observing 
methods, the data processing, and the amplitude calibration is 
described in Kovalev et al. ( 1999). 

All single-dish spectra obtained with the Effelsberg 100-m 
and RATAN-600 radio telescopes are presented in Figure |2l 
Here, no indication of a flare or strong difference/variability 
between August and September can be noted. However, the 
December 2008 spectrum shows the beginning of a bright ra- 
dio flare with a clear spectral steepening towards higher fre- 
quencies {1/ > 10 GHz). This is in good agreement with the 
strong flux density increase seen in the 15 GHz light curve 
during November/December (Fig|7]). 

Detailed radio images at sub-milliarcsecond scale of the 
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Fig. 6. — Long term radio flux light curve at 15 GHz obtained by the Owens 
Valley Radio Observatory (OVRO) 40m dish radio telescope (filled circles), 
showing the rising part of a radio outburst started in late November 2008, 
i.e. almost 4 months after the gamma-ray outburst detected by Fermi. The 
fill diamonds represent the flux measurements performed by the Metsahovi 
14-m radio telescope at 37 GHz (right y-axis scale), confirming the start of a 
radio outburst at a higher frequency. The scaled LAT daily light curve on the 
same period is reported for comparison (small light grey bars). 
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Fig. 7. — Variable broad-band radio spectra observed with the Effelsberg 
100-m and RATAN-600 radio telescopes simultaneous to the LAT data. His- 
torical RATAN-600 data (grey open triangles) and archival data from the lit- 
erature until March 2008 (grey open circles) are shown in the background for 
comparison. 

PKS 1502-1-106 superluminal jet were obtained during three 
epochs in 2008 with Fermi ah^eady in orbit: on June 25, Au- 
gust 06 (during the maximum peak of the 7-ray outburst), and 
November 19. These observations were performed as part of 
the MOJAVE monitoring program conducted with th e Very 
Long Baseline Array (VL B A) at A = 2 cm (Lister & HomanI 
120051; iLister et ani2009al) and provided useful high resolution 
total intensity and linear polarization images. These VLBA 
images close to the 7-ray flare are reported and compared to 
the map obtained one year earlier in Figure H] The highest 
integrated flux density val ue since the beginning o f the 2 cm 
VLBA monitoring in 1997 (iKeUermann et al.l2004 ') was mea- 
sured on November 19, 2008 (DOY 324) as F15 ghz =2.0 Jy, 
with a peak intensity of 1.6 Jy beam~^. These values are sig- 
nificantly hig her than the typical l evel of 1.3 Jy repo rted in 
the program (iKovalev et all 120051: iLister et al.1 l2009allbl) and 
indicate a radio flare happening in the source VLBI core. 
The core flux density and brightness temperature raised to 
higher values as well which means that the flare happens in 
the VLBI core, as expected. These finding are in good agree- 
ment with the single-dish results presented above (Figure |6] 
and Figure |7ll. The second relevant feature is the direction 
of the electric vector position angle (EVPA) in the core re- 
gion, which rotated between the 2007 and 2008 epochs by 
90 degrees most probably indicating an opacity change — a 
precursor of an outburst in the VLBI core. 

In summary, our single-dish and VLBI radio monitoring 
of PKS 1502-1-106 simultaneous to the Fermi LAT observa- 
tions has revealed (i) no significant radio 15 GHz variabil- 
ity during the strong LAT 7-ray flare seen in August 2008, 
and (ii) a strong radio flare which becomes clearly visible at 

> 10 GHz during October/November 2008 (Figure |6]l with 
a rise phase lasting for at least 20 days. If this flaring behav- 
ior is associated with the bright 7-ray flare of August 2008, 
a delay of more than three months (^ 98 days if the lag 
between the starting days of the 7-ray outburst, DoY 218.2, 
and the 15 GHz outburst, DoY 316.7, is considered) could be 
expl ained by opacity e ffects in the core region of the source 
(e.g.. AUer et al.|[r999h . Although we cannot exclude the pos- 
sibility of radio activity at 15 GHz during the OVRO 40-m ob- 
servations outage (September 3 - Otober 23 2008) that might 
also be associated with the gamma-ray flare of August 2008, 
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Fig. 8. — Total intensity and linear polarization images observed by the 
VLBA at 15 GHz as part of the large Fermf'-supporting MOJAVE program. 
Naturally-weighted total intensity images are shown by black contours. The 
contours ai'e in successive powers of two times the base contour level of 
1.0 mJy beam^^. Electric polarization vector directions are indicated on the 
left hand side by blue sticks, with their length being proportional to the polar- 
ized intensity. Linear fractional polarization is shown on the right hand side 
overlaid according to the color wedge. 



OVRO data before and after the outage as well as September 
and October Effelsberg 100-m data are consistent with very 
little change in the 15 GHz flux over this time period. The 
three flux measures at 37 GHz obtained at the Metsahovi ra- 
dio observatory during the OVRO outage, confirms this trend 
with very little variability during this period. However, the 
observed radio flare could also be associated with, e.g., the 
more recent, prominent variability seen in the LAT 7-ray data 
during November/December (Figured DoY~' 320-333). A 
more detailed analysis of such possible correlations and the 
source's overall radio/7-ray behavior seen with LAT and si- 
multaneous radio observations over a longer period of time 
will be the subject of a subsequent work. 

4. ARCHIVAL MULTIWAVEBAND DATA 
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A full multiwavelength an alysis dedicated o n PKS 
1502+106 is available only from lGeorge et alj (119941) . where 
old archival and broadband radio to X-ray data obtained in 
1993-94 were presented. Data and analysis on PKS 1502+106 
reported in other papers are mostly limited to the radio regime. 
In order to compare our multifrequency findings with the past, 
and to form a more complete characterization of this blazar, 
we briefly present here, for the first time, results from un- 
published past observations by INTEGRAL, XMM-Newton, 
Swift, and Spitzer space telescopes performed in 2001, 2005 
and 2006. 

4.1. INTEGRAL observations in 2006 

The sky region containing PKS 1502+106 was observed in 
2006 by IBIS (83ksec, MJDs 53760.4 to 53762.4, Jan.25-27, 
2006), and a new soft 7-ray source (IGR J 15039+ 1022) was 
detected with a flux density of 1.6 mCrab in the 18-60 keV 
energy range (corresponding to 1.2 x 10~^^ erg cm^^ s^^, 
see ATel #1652). This IBIS source was identified with Mkn 
841, a Seyfert galaxy known to display a well detected high 
energy cut-off around 100 keV (Petrucci et al. 2002), making 
it unlikely to emit in the 7-ray domain. The angular distance 
of PKS 1502+106 from IGR J15039+1022, - 11', points to 
a clear non-detection during this Jan. 2006 INTEGRAL ob- 
servation, while a 2a upper limit for PKS 1502+106 of 0.7 
mCrab in the range 18-60 keV (0.52 x 10"" erg cm"^ s"^) 
is inferred. 

4.2. XMM-Newton and Swift 

Four serendipitous, unpublished XMM-Newton X-ray ob- 
servations of PKS 1502+106 by the EPIC (MOS detector 
only) camera are available as the source was in the frame of 
the target Seyfert galaxy Mkn 841. PKS 1502+106 was al- 
ways on the border of the MOS chips and out of the PN frame, 
and therefore subject to low X-ray statistics, regardless of its 
intrinsic brightness. The four X-ray EPIC-MOS observations 
(three in 2001 and one performed in Jul. 17, 2005, see Table|3] 
for an analysis summary) do not show variations in the 0.2-10 
keV photon index, while the 0.2-10 keV flux intensity var- 
ied by a factor of a few (in the range 3.5 — 6.8 x 10" ^"^ erg 
cm~^ s~^in 2001, comparable to the lower states observed by 
ASCA, and a mildly active state with 1.1 x 10"^^ erg cm"^ 
s-Hn July 2005, TableO. 

PKS 1502+106 was also observed twice in the past with 
the Swift-XKI as a fill-in target (TargetlD: 36388), showing a 
fainter X-ray flux (0.02 counts/sec) than the flux recorded in 
the August 2008 campaign observations (0.05 counts/sec). 

4.3. Spitzer observations and the multifrequency behavior on 
Jul-Aug.2005 

In the past, only upper limits in the far-/near-infrared 
bands by IRAS were available for PKS 1502+106 
(lNeugebaueretal.1 119861 and Figures ^ and [Ip. PKS 
1502+106 was observed serendipitously in the mid-infrared 
band for the first and only time by Spitzer on August 13, 2005, 
09:10-09:18 UT (PID 117, A PR Reques t Key 5011456). The 
Infrared Spectrograph (IRS, Houck et al . 2004), low resolu- 
tion (R=60-130) module, recorded the mid-IR spectrum from 
5-14 /im (shown in Fig |9] with the optical SDSS spectrum 
and near-IR photometric data point). High resolution IRS 
module spectra were also taken, but could not be used since 
there were no accompanying background observations. The 
Short-Low (SL) coadded 2D spectra were reduced using 
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Fig. 9. — The unique observed mid-IR spectrum in the range 5-14 /xm 
obtained by the Spitzer Infrared Spectrograph (IRS) low resolution {R = 
60 — 130) module, in August 13, 2005. The position of redshifted Brackett 
emission lines and PAH hne are indicated, even if they are not detected in 
the IRS spectrum that is consistent with a simple power law model. In addi- 
tion the optical spectrum by the SDSS on April 23, 2006, and high precision 
JHKg photometric flux measurements of the Kitt Peak National Observatory 
(KPNO) 2.1m telescope of March 29, 2002 are also reported. 

the standard Spitzer IRS pipeline (ver S17.2). Background 
was subtracted using the two nod positions along the slit. 
The spectra were extracted and flux-calibrated with SPICE 
ver. 2.1.2, in a standard, expanding point-source aperture. 
The two spectral orders match well at 7 /im, indicating a 
well-pointed observation. The mid-IR continuum of PKS 
1502+106 rises to the near infrared and appears to be rather 
featureless, consistent with pure synchrotron emission (a 
power law F^, cx v^'^-^ in the 5-14 ^m range, Fig|9l). The 
wave-like deviation of the data can be simply explained as 
wavelength-dependent spectrograph slit losses, while the 
red-shifted Brackett emission line series (like the 3.3 /im 
PAH feature) falling in this wavelength range are indicated, 
but they are not well detected in the spectrum. This IRS 
spectrum is similar to other blazars which have been observed 
by Sp itzer, including BL L ac and 3C 454.3 dLeipski et alj 
2008; Ode et al. et"al]|2009l) . while flie near IR (J,H,K) flux 
data from the Kitt Peak National Observatory (KPNO) 2.1m 
telescope Watanabe et al. (2004) reported in the same figure 
indicate a lower flux state and steeper H-Ks (1.65-2.15 /^m), 

TABLE 3 

Analysis summary of the EPIC-MOS instrument observations 
(Jan. 2001 and Jul. 2005) on board of XMM-Newton. 



Obs. id. (date) Best-fit Model and Parameters 

ObslD 0070740101 Power Law 

(Jan 13, 2001, 09:20 UTC) rx=1.6±0.2 
X?=l -69/11 

Fo.2-10 keV = 3.5 X IQ-" erg cm~^ s'^ 

ObsID 0070740301 Power Law 

(Jan 14, 2001, 00:30 UTC) rx=1.7±0.2 
X?=l-27/17 

Fo.2-10 keV = 6.8 X 10-13 erg cm"^ s-^ 

ObsID 1 1 29 1 020 1 Power Law 

(Jan 13, 2001, 04:58 UTC) rx=1.6±0.2 
X^=105/8 

Fo.2-10 fceV = 3.6 X 10~" erg cm~^ s'^ 

ObsID 0205340401 Power Law 

(Jul 17, 2005, 06:32 UTC) rx= 1.69 ± 0.08 
X^=0.99/58 

Fo.2-10 keV = 11-0 X IQ-^^ erg cm~^ s~ 
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spectral index {a(^H-Ks) = 1-66). 

The July-August 2005 SED assembled with these XMM- 
Newton (July 17, 2005) and Spitzer (August 13, 2005) ob- 
servations (joined with a couple of radio-optical observations 
in these two months by the Metsahovi, RATAN and Catalina 
observatories) is consistent with a low or mildly active stage 
that can be explained by a simple SSC model (inset plot in 
FigHD- Based on these data there are therefore no hints for 
a deviation from a SSC scenario in this flat spectrum radio 
quasar. In any case, this SED cannot provide significant con- 
straints, because it contains observations obtained over about 
two months and does not include any 7-ray data, aside from 
the older (19 92) EGRET upper limi t (calculated with the pre- 
scriptions in [Thompson et al.|[T996l) . 

XMM-Newton observations pointed out an (0.2-10 keV) 
photon index that shown almost no variations, with val- 
ues (Fx ^ 1.7) in agreement with the previous X- ray ob- 
serva ti ons performed by RO SAT and ASCA (Geor geet alj 
119941; IWatanabe etall 120041) . The integrated X-ray flux 
Fo.2-10 keV = 1.0 X 10~^^ suggests a mildly active state dur- 
ing this observation, comparable to the Swift XRT spectrum 
observed in the days soon after the LAT outburst (red filled 
symbols, data points on FigfTTTl. More complicated and pos- 
sibly more accurate emission models, beyond the SSC, can 
be therefore investigated for the first time only thanks to the 
Fermi multifrequency campaign of August 2008, whose SED 
is described in the next section. The Swift XRT simultane- 
ous spectrum of Aug. 2008 has a slightly harder spectrum with 
respect to these archival XMM-Newton observations while 
the XRT flux was about one order of magnitude greater than 
the flux observed in the past by XMM-Newton, ASCA and 
ROSAT. 

5. DISCUSSION 

5.1. Gamma-ray outburst and longer-term variability 

During the first several months of the LAT survey, PKS 
1502+106 was one of the brightest, as well as the most lu- 
minous, blazar in the MeV-GeV band. The threefold flux in- 
crease in < 12 hours between Aug. 05 and Aug. 06, 2008 
(DoY 218-219, Fig. |2]l constrains the rest-frame size (R') of 
the flai-ing region: R' < cAtV/{l + z) ~ 6.8 x 10^^ cm 
(where V = l/(r(l — f3cos9)) is the macroscopic Doppler 
factor, r the bulk Lorentz factor and 6* ~ l/F the angle of 
sight). The value 2? ~ F ~ 0a„ r> = 14.8 is assurn ed from 
MOJAVE VLBA measurements dLister et alj|2009bl) . where 
Papp is the kinematic apparent jet speed in units of c. In this 
high speed regime an upper limit on the viewing angle can be 
also estimated: 9 < 2 arctan(l//3Qpp) = 7°. 7. This is con- 
sistent with an independent estimation of the Doppler factor 
based on the 37 GHz flux variability, made in iHovatta et alj 
(120091) : a factor 2?„t„. = 12 and Papp.var = 14.6, with bright- 
ness temperature T^r — 8.7 x lO^'^ °K, and angle of sight 
e = 4.7°. 

Assuming the "concordance" cosmology (Sect. [T]and z = 
1.839), the luminosity distance of PKS 1502+106 is d-L = 
14.2 Gpc, and the inferred, apparent and isotropic, monochro- 
matic luminosity at Eq = 100 MeV of PKS 1502+106 during 
the outburst phase is i_B>ioo — 47rci| • (F — 1) EoFe>Eo — 
1.1 X 10"'^ erg s^^,, where the average flux Fe> Eq = 2.91 x 
10-6phcm-2 s-i in the outburst interval (DoY 2008: 218.95 
- 224.0) is used. The bolometric luminosity is expected to be 
even higher than this value, since the measured LAT spectrum 
appears to be beyond the peak of the high energy component. 



and therefore this LAT blazar has probably one of the highest 
Ly/At ratios (2.5 x 10"*^ erg s^^) known in the MeV-GeV 
regime. 

Relativistic motion provides a solution for the problem of 
intrinsic excess absorption by pair-production in powerful 
7-ray sourc es like PKS 15 02+106 (see, e.g. iMat tox et al.l 
1993; Mad eiski et al.lfT996l) which have a significant Lblr- 
Adopting the flux tripling time scale of the outburst rise (i.e. 
A< = 12 hours), and the outburst state averaged X-ray flux 
(Fo.3-10 keV = 2.18to i2 X 10-12^ Table|2]l at the observed 
photon frequency vx = 10^^ Hz, the minimum Doppler 
factor T) required for the photon-photon annihilation optical 
depth to be t^^ < 1 can be estimated. Using the derived 
relation 1 = r..^^ ~ ard^Fx / {iAt ExV^) and taking 
the region size R — cAtV/ (1 + z), the source-frame photon 
energy E'^ = (1 + z)hvx/'D and the intrinsic X-ray lumi- 
nosity L'x Andl'D^'^Fx we obtain V > 7.7 (omitting the 
requirement of co-spatiality of the X-ray and 7-ray emission 
regions relaxes this limit). This is in agreement with the val- 
ues found from radio flux-structure observations and with the 
SED modeling parameters found. 

The asymmetry of the August 5-6 7-ray outburst can sug- 
gest a more complex emission geometry than a simple one- 
zone model. The temporal structure — ^ 0.5 day rise, fol- 
lowed by a ^ 4.5 day decay where a ^ 2.5 day intermediate- 
level plateau is likely observed — implies particle accelera- 
tion and cooling times that are greater than the light cross- 
ing time, i.e., tinj,icooi > R/c, where these quantities are 
evaluated in the jet comoving frame. A synchrotron self- 
Compton (SSC) emitting blob in the jet should be relatively 
confined (< 0.01 pc), although relativistic beaming would 
permit the region to be as much as an order of magnitude 
larger The hinted intermediate plateau could mark a twofold 
active region, and two SSC emitting components. Descrip- 
tions making use of a multi-zone SSC or multi-emission 
component SSC (second order, superquadratic components) 
are reported, for example, i n Sokolov & Marscher (2005); 
lGeorganopoulosetal.1 (l2006h ; iGraffetal.. (.2008) . On the 
other hand, if the injection of relativistic electrons is impulsive 
and repeated, some single-zone SSC models already predict 
platea us during an outburst (see^e.g., Chiaberge & Ghisellin^ 
119991; iBottcher& Chiang! 120021: ISokolov et al.l 120041) . At 
lower frequencies (IR-optical-UV), where cooling times are 
longer, the electron distributions corresponding to different 
injections can build up, and the memory of the individual in- 
jection phases can be lost, providing a smoother decay with- 
out intermediate plateaus (as shown in X-ray, optical-UV light 
curves by Swift, Fig. |4|i. Apparent delays (like the 4-day lag 
hinted by UVOT data) can also be explained within this sce- 
nario. 

On the other hand this asymmetric (fast rise, slower decay) 
shape of the 7-ray outburst can be also an evidence for a dom- 
inant contribution by Comptonization of photons produced 
outside the jet (Sect. 15.31 and 15.4b during this event, as pre- 
dicted for example in lSikora et al.. (.2001 ). Gamma-ray flares 
produced by short-lasting energetic electron injections and at 
larger jet opening angles are predicted to be more asymmet- 
ric showing much faster increase than decay, the latter deter- 
mined by the light travel time effects. 

The "harder when brighter" trend of the gamma-ray pho- 
ton index during the outburst (right inset panel of Fig. |2|i, 
hints a narrow hysteresis evolution of the spectral index 
against the flux, a signature produced by non-thermal cool- 
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ing and high to low energy p r opagation of the electron in- 
iection rate dKirk et alj [l998t iGeorganopoulos & Marscheii 
Il998l: iBottcher & Chiang|l2002D . The photon index extracted 
with a power-law model over the daily bins, was found quite 
scattered irrespective of the flux level in the remain part of the 
light curve following the outburst. 

The outburst of Aug. 2008 appeared to have ignited an en- 
during relaxing state of 7-ray brightness and activity, dur- 
ing the following four months (Sections 12.31 and 15.11 1. The 
1//" (with a ^ 1.3) PSD points out a general fluctua- 
tion mode placed between a pink-noise (flickering) and a 
random-walk, (Brownian motion or brown noise), staggered 
by some rapid flares, similar to the long-term variability of 
blazars observed in r adio and optical \yavebands (for example, 
Hughes etal."1992'; 'Ciaramella e t all l2004t iTerasranta etll] 
2005; Ciprinietal. 2007b; Hova tta et al.1 boOTh . In con- 
trast, this variability mode is rather different than the full 
Brownian regime shown by the short-term (intra-hour reso- 
lution) hght curve of di e large TeV outburst of PKS 2155-304 
jAharonian et alj|2007h . This PSD indicates the occurrence 
frequency of a specific variation is inversely proportional to 
its strength, as found in processes driven by stochastic relax- 
ation, and rapid flares/outbursts (phenomenologically called 
intermittence), are not occasional events produced by physical 
processes of different nature with respect to the mechanism 
responsible for the long-term flickering, but can be considered 
as statistical tails of the same dynamic process, possibly con- 
nected to disk or jet instabilities, to viscosity and magnetic 
turbulence, or to inhomogeneities and shocks (for example 
iBegelman et al. 1984 ). Even in the case of quite nonhomoge- 
neous structures, the jet is seen under a small viewing angle 
{9 < 7°. 7 here), therefore flaring events and variability trends 
are the result of emission components originating from differ- 
ent regions, excluding, in most cases, causality. 

5.2. Gamma-ray and radio connection 

The absence of significant radio 15 GHz and 37 GHz vari- 
ability before, during and immediately after the LAT outburst 
(Fig. |6]l, and strictly correlated to the near-lR to 7-ray out- 
burst, has consequences for opacity at lower radio frequen- 
cies. However, the beginning of the strong radio flare seen in 
the Oct.-Dec. 2008 data of Fig. |6]and Fig. ^\ could be as- 
sociated with the period of increased 7-ray activity seen after 
-DoY 330 (Nov. 25, 2008, Fig. Ell. If the beginning of the ra- 
dio outburst in Fig. |6]is associated with the LAT outburs t, thi s 
would imply a delay of more than about 98 days (Sect. 13.3b . 
meaning that the radio emission region becomes transparent 
at the cm and longer wavelengths significantly after the out- 
burst seen at wavelengths where the source is truly optically 
thin (mm, IR, optical bands). 

The rotation direction of the electric vector position angle 
(EVPA) in the milHarcsecond-scale core of PKS 1502H-106 
between 2007 and 2008 was a possible precursor signature 
of an outburst that occurred in the core (Fig|8]l. The elec- 
tric polarization vectors appear well aligned to jet axis in the 
June 25, 2008 map and even more in the Aug. 06, 2008 map. 
This latter map represents an unprecedented example of radio- 
structure snapshot truly simultaneous to the peak of a MeV- 
GeV outburst (Fig|8]l. About 3 months after the outburst, on 
Nov. 11, 2008, the alignment is again decreasing. These in- 
teresting findings are in agreement with the scenario that as- 
sumes very bright 7-ray flares occurring after the ejections of 
superluminal radio knots, with accompanying increases in po- 
larized radio flux, and a field ordering and alignment with re- 



spect to the jet axis. A correspondence between variations in 
polarization direction and intensity in different bands at par- 
sec scales can help localize the primary site of the high energy 
emission. We found that the prominent 7-ray outburst and 
the possible 3-month delayed radio-outburst are likely pro- 
duced in the core of the parsec-scale jet of PKS 1502H-106. 
More details on connections between VLBl radio structures 
and high energy emission, and their interpreta tion are reported 
recently in Jorstad et al. (2007); d'Arcange lo et al.1 (12007); 
'Marsc her et' al.i (i2008i) ; iKovalev et ajj (i2009i) . To test if this 
is indeed the case for the discussed gamma-ray flare, VLBA 
monitoring of 1502H-106 is continued at 15 and 43 GHz. 

5.3. Gamma-ray, X-ray and UV-optical cross correlations 

During the EGRET era a similar degree of simultaneous 
X-ray and 7-ray m onitoring was achieved, for instance, fo r 
the blazar 3C 279 ( Wehrle et al.|[T998^. 'Hartman et aLll^OOT), 
a typical 7-ray powerful FSRQ resembling PKS 1502H-106. 
The correlation found in PKS 1502H-106 between the 7-ray 
flux and the X-ray and UV-optical emission during the 16- 
days of Swift follow up, is evidenced by the multi-panel and 
multifrequency light curves in Fig. HI The apparent linear cor- 
relation between the emission in the LAT band and the X- 
rays (Fig.[Tol left panel) suggests that the observed X-ray and 
MeV-GeV photons may be part of a single SSC component 
(continuous line, labeled "SSC2", on FigfTTI). This could also 
explain the nearly constant X-ray (0.3-10 keV) spectral in- 
dex observed during the outburst with respect to the post-flare 
phase. 

The 7-ray peak power observed during the outburst was 
vFiy ^ 10"*^ erg s^ , and decreased to less than about 2 x 10^^ 
erg s^^ in the post-flare period, meaning a difference of at 
least a factor of ^ 6 (Fig. fTTT i. The difference between these 
two emission states in X-ray output was instead a factor of 
^ 2.5, and a factor of ^ 3.5 in the optical-UV band (frequen- 
cies above the synchrotron peak). These sub-quadratic dif- 
ferences between the outburst and the subsequent state could 
be still described by SSC descriptions. The strict correlation 
between the X-ray flare and the 7-ray outburst supports, at 
least, the dominance in the SED (FigfTTli of the SSC, in-jet, 
emission fr om radio to X -ray bands in agreement with results 
detailed in Sikora et al.l (|2001). A simple single-zone SSC 
model has nevertheless problems explaining and reproducing 
the large 7-ray dominance observed during the outburst, and 
could require very sub-equipartition magnetic fields. 

Fig. E] and the right panel of Fig. [TO] show the flux mon- 
itored by UVOT (U-band reported there, but similar results 
are found in the other 5 filters), indicating the same strong 
correlation with respect to 7-rays, but with a possible time 
lag of about 4 days (also hinted at by the peak on the dis- 
crete cross-correlation function, DCCF, inset plot). This pos- 
sible time lag can be reasonable only by assuming that the 
optical-UV brightness during the start and rise of the gamma- 
ray outburst (i.e. during DoY 216-219) was comparable or 
lower than the flux observed during the first UVOT observa- 
tion (performed on DoY 220.8). Relative delays between the 
synchrotron (our Swift UVOT data) and the inverse Compton 
(our Swift XRT and Fermi LAT data) counterparts, are dom- 
inated by energy stratification and geometry of the emitting 
region. In SSC-dominated zones the synchrotron emission 
is co-located with the inverse-Compton production site, and 
such 4-day optical-UV to 7-ray lag would depend on light- 
travel time effects in the the emitting region (characterized 
by size R and viewing angle 9), and by the particle cooHng 
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Fig. 10. — Left panel: the gamma-ray flux measured by Fermi LAT versu.s 
the X-ray flux measured by Swift-'KKT. The cross-coiTelation, without lags 
is well displayed. Right panel: the LAT gamma-ray flux versus the UV flux 
measured by Swift-\]\OT in the representative U-band, and the discrete cross 
correlation function diagram (DCCF, inset plot) between the gamma-ray flux 
and the U-band flux. Here a clear correlation with a 4 day delay of the UV 
emission flare with respect to the gamma-ray flare is suggested. 

time and decay time of the synchrotron and inverse Comp- 
ton (IC) emissions. When the region is not so compact, both 
decay times can be comparable to the apparent light cross- 
ing time R/c, and significant shifts between light curves at 
these different energy bands are expected. Synchrotron flares 
(our optical-UV data) can be delayed with respect to the IC 
7-ray flares (e.g. lChiaberge & Ghisellini|[l999l:ISokolov et al] 
12004). Rather long time lags such as this could also be ex- 
plained by a prolonged disturbance traveling down the jet. 
The disturbance, being radially inhomogeneous in both den- 
sity and velocity, could induce shocks and collisions leading 
to the formation of two adjacent emission zones with similar 
properties (multi-zone SSC), thus explaining the flare shape 
asymmetry and the intermediate-level plateau during the de- 
cay (Section lOT l. 

Unlike SSC emission, outbursts dominated by inverse 
Compton external radiation involve a constant field of seed 
photons, and are not delayed by light-travel time of pho- 
tons. The resulting frequency stratification behind the shock 
front could extinguish an outburst first at the highest ener- 
gies, then progressing to low er frequencies as time advances 
(ISokolov & Marscheill2005l) . and energy-dependent radiative 
losses induce delays in the declining part of the emission pro- 
duced by the lower energy particles. The variability patterns 
shown by our data could rule out radiative cooling as the only 
mechanism for causing the delay, since the rising part of the 
optical emission follows the rising part of the 7-rays. In ad- 
dition, the decay time scale in the UVOT data seems com- 
parable to that seen in the LAT, and not longer as required 
by energy-dependent dominated cooling. A reasonable mix- 
ture of SSC and an extra-contribution by Comptonization of 
the BLR photon field (external jet origin) by the same popula- 
tion of energetic electrons (ERC, e.g. Sikora et al. 1994, 20011 
iDermer & SchUckeiser 2002: GhiseUini & Tavecchio 2009aR 
could be more adequate to explain also the SED, and such 
SSCh-ERC scenario is still compatible with the X-ray to 7-ray 
con-elation found in PKS 1502H-106. 

As seen in Section 13.21 the first of the fast, minor 7-ray 
flai-es occurring in Sept. 2008 (DoY 248-251) after the big 
outburst, was well observed by the Kanata-TRJSPEC tele- 
scope in both J and V bands, with no significant time lag. 
The lower intensity and duration of the flare, the limited tem- 
poral resolution of the data, and a possible dominance of the 
SSC process during this episode can explain this difference 



with respect to August's big outburst. The match between the 
MeV-GeV and optical-near-lR flare, was also crucial for the 
firm identification with PKS 1502H-106 (Sect. |Z2T |, and was a 
reasonable confirmation of the source-intrinsic nature of this 
variation seen in the LAT light curve (Fig.[2]and Fig.|4]i. These 
V and J flux measurements, obtained after the conclusion of 
the outburst phase but on a period longer than the Swift moni- 
toring (Fig.|4ll, appear well correlated, though the V — J color 
was less pronounced than the multi-band near-lR colors re- 
ported in'Watan abe et al.l (12004 ) and Fig. |9] while the degree 
of optical polarization remained almost constant. 

5.4. Spectral energy distribution 

Variability is a powerful diagnostic to investigate blazar 
physics, but represents also a supplementary problem for the 
analysis of broad band spectral energy distributions (SED), 
where model constraints are provided by simultaneous and 
well time-resolved multifrequency data. In the previous sec- 
tion we indicated that dominant synchrotron and SSC (in- 
jet) mechanisms can explain the radio-to-X-ray emission and 
correlations, whereas the origin of the high-power MeV- 
GeV bolometric emission could be better constrained by ex- 
ternal Comptonization of the radiation from the broad line 
region (BLR), as invoked by (Ghisellini et al. 2009b), and 
found already in similar FSRQs of the EGRET era (e.g., 
ISokolov &Marscheril2005L Sikora etaP 12008. 2009) . The 
featureless mid-IR continuum observed by Spitzer-IRS in 
2005 (Fig. |9]l, supports the hypothesis of a prominent syn- 
chrotron emission by the jet, controlling the lower energy 
component. It could also reflect the lack of detectable am- 
bient dust radiation, thereby supporting the idea that the ERC 
dissipation occurred within the BLR for this LAT o utburst, in 
agreement with prescriptions of (ISikora et al.ll2002h . 

PKS 1502-1-106 might be considered a blazar peaked at the 
border of the MeV and GeV bands (a peak around 0.4-0.5 
GeV is suggested by the curved model fit of the 4-month spec- 
trum, Section|231i. In other words, PKS 1502H-106 is Ukely 
at the border of the family of BLR-dissipated FSRQs and 
circum-nuclear ambient/torus dust-dissipated FSRQs, with an 
important SSC power output from radio-to-X-ray bands (Sec- 
tion 15. 3t , as depicted by the multiband correlations, the ab- 
sence of hints for a bulk Compton feature produced by cold, 
adiabatically cooling electrons (as obser ved in PKS 1510-089, 
iKataoka et ani2008t lAbdo et al.ll2009fl) . flie lack of evidence 
for a blue bump. 

The medium or high black hole mass of this blazar 
(likely in the range 0.5 — 1 x 1 0^ Mq, as calculated by 
D'Elia et al. 2003; Li u et al.l 120061 respectively) could con- 
note a BLR radiation field and inverse Compton dissipation 
moderately stronger than the magnetic energy density and the 
SSC luminosity in the gamma-rays, as reporte d in Fig. [TTI 
The estimated accretion rate is 2 Mq yr~^ (iD'Elia et al.l 
l2003h . In this case we have Lerc / Lsyn = t^Bijj/C^s — 
LblrP^ I (47rci?^^^C/s) where Ugj^j^ and U b are the elec- 
tromagnetic energy density in the BLR and the magnetic 
field energy density in the jet blob, respectively. With a 10^ 
Mq mass the Schwarzschild radius of the supermassive black 
hole is Rs = 2GA//c2 ~ 3 x 10^^ cm. The bolomet- 
ric luminosity of the BLR in PKS 1502-1-106, evinced by the 
Mgll emission line profile, is Lblr ~ 3.7 x 10*^ erg s^^ 
dLiu et al.ll2006l) . implying that the BLR is located at a radius 

Rblr = lO^VVlO"*^ - 2 X 10^^ cm. We can expect rela- 
tively smaller magnetic fields in the dissipation region and an 
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Fig. 1 1 . — Major panel: Overall radio-to-gamma-ray spectral energy distribution (SED) of PKS 1502+106 assembled with data from the 2-week Fermi LAT 
multifrequency campaign of August 2008. The two representative, time-averaged, states for this campaign (time intervals outHned in Fig|2), i.e. the high state 
(Aug. 05-10, 2008, DoY 2008: 218.95 - 224.0; blue filled circle symbol), and the post-flare (intermediate brightness) state (Aug. 1 1-22, 2008; DoY 2008: 224.0 - 
235.42; red filled square symbol), are represented along with their SSC and ERC model attempts. Ai'chival non-simultaneous data (including the Jul. -Aug. 2005 
data, and the whole 4-month cumulated LAT spectrum, open orange circle symbols and strip) are reported in the background for comparison. The high Compton 
dominance and gamma-ray bolometric luminosity reached during the outburst is evident. The SED of the outburst state is reported with a supeiposed one-zone 
"SSC+ERC" model fit attempt (where the SSC dominates the radio-to-X-ray SED and the ERC produces the gamma-ray component, blue dashed line), while 
the SED of the post-flare state is reported with two possible models superposed: the same SSC+ERC modeling (SSC radio-to-X-ray and ERC for gamma-ray 
band, red dotted line) and also a pure one-zone SCC attempt (labeled with "SSC2", a second type of stand-alone SSC model for the entire radio-to-gamma-ray 
SED, continuous red line). Inset panel: The July-August 2005 non-simultaneous SED (filled orange circles) assembled with the XMM-Newton EPIC-MOS data 
from Jul. 17, 2005, the Spitzer IRS observation from Aug. 13, 2005, radio flux data from Metsahovi and RATAN radio observatories, and optical data from the 
Catalina Sky Survey. Archival non-simultaneous data (open circles) are reported in background. This two-month averaged SED is consistent with a low or mildly 
active stage (suggested by the X-ray flux above 1 X 10~^^ erg cm~^ s~^) and can still be explained by a pure one-zone SSC model. 



ERC, external-jet process that is the dominant cooling mech- 
anism when the source is in a high 7-ray state, as suggested in 
Section ISTl In fact a SSC-only model would require a mag- 
netic field in very sub-equipartition conditions, even if only 
the BLR fraction lying within the varying region's beaming 
cone could contribute to gamma-rays. The simultaneous SED 
corresponding to the outburst state reported in Fig. [TT] cor- 
roborates the possible needs for an ERC contribution, show- 
ing a large {L^/Lopt ^ 100) gamma-ray dominance over the 
synchrotron component. The SSC-i-ERC model can be con- 
sidered in some way in agreement with the redshift value of 
this blazar, if we assume the FSRQs BL Lacs cosmolog- 
ical evolutionary scenario with dimming jet power, possibly 
related to star formation rate or the the far-IR/submillimeter 
luminosity density (Dermer 2007). 

In Fig. [TT] we report the data of the two SED states (the 
"outburst state" of Aug. 05-10, 2008 (blue points) and the 
"post-flare" state (red points), characterized by a intermediate 
luminosity, of Aug. 11-22, 2008 (both intervals are outlined 
in Fig. |2]i, in conjunction with the unpublished archival and 
literature data collected for comparison. Above the two si- 
multaneous SEDs of Aug. 2008, SSC and ERC "strawman" 
models are reported for each state (SSC-i-ERC modeling, blue 
lines for the outburst period, and SSC-i-ERC plus a SSC-stand- 
alone modeling, red lines, both for the post-flare period). 



The pure SSC stand-alone model for the post-flare state (in- 
dicated in Fig.[TT|as a version 2 of the first order SSC model- 
ing with the label "SSC2", continuous red line, spanning from 
radio to gamma-ray bands) implements the temporal evolu- 
tion of the synchrotron and SSC spectral components in a sin- 
gle flaring blob within the jet, where a population of acceler- 
ated electrons having a power-law with exponential cutoff dis- 
tribution (dN/dE oc £;-Pe~(-^/-^""")), is instantly injected 
(for more det ails on the an alytical and numerical model see, 
for example, ICiprinil 120081) . This form for dN/dE is plausi- 
ble in prese nce of time-dependent accelera tion or radiative- 
loss limits jWebb et alj [19841 lDrurvl[T99l . The post-flare 
(lower) state is reproduced by setting the electron energy in- 
dex p = 1.77, minimal and maximal electron Lorentz fac- 
tors 7min = 200, 7ma£c = 8 X 10^, compactncss injected 
iinj = LinjCTT / {RnieC^) = 10^^, radius of the emitting 
region in the comoving frame R = 6.5 x 10^^ cm, a bulk 
Doppler factor 2? = 8 and magnetic field intensity B = 0.024 
G. If the X-ray spectrum is produced by electrons that cool on 
timescales longer than the light crossing time, R/c, the X-ray 
spectral index would be ax = (p— l)/2~0.4, a value quite 
similar to the averaged value {ax = 0.45±0.03) measured by 
Swift XRT for the post-flare interval (bottom segment of Table 
[2]). The rather large size of the emitting region (R ~ 0.2pc) 
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and the very sub-equipartition magnetic field and reduced V 
suggested by this "SSC2" (stand-alone and first order radio- 
to-gamma-ray SSC) model attempt indicate that even for this 
lower brightness state, an ERC contribution can be reason- 
able. Furthermore, the value of the optical-UV spectral in- 
dex for this post-flare state {auvoT — 1.9 ± .3) is softer 
than the averaged index (3/2) expected jChiang & Bottcher 
120021) . when a totally SSC-dominated loss is considered. The 
huge Compton dominance of the outburst state might preclude 
a first order SSC-stand-alone model attempt for such SED 
state. On the other hand for high-energy 7-ray loud blazars 
like PKS 1502+106 superquadratic variations produced by 
higher ord er scatters in SSC are predicted by oth er models 
(iGeorgano poulos et all 120061: iPerlman et all 120081) and Sec- 
tion 15 . 1 1 in fact in luminous blazars, the optical depth to 
Compton scattering (related to particle density) increases and 
second and higher-order scatters become more important. In 
more Compton-dominated objects these scattering reactions 
can dominate the energy output from the SSC process, and 
can produce superquadratic behaviors during big flares. 

As mentioned above a "SSC plus ERC" model descrip- 
tion, joining a SSC radio-to-X-ray component and a ERC 
gamma-ray component, to both the post-flare and outburst 
states (FigHT] dashed blue and red dotted lines) is plausible 
as well for these SEDs. The description of an initial version 
of this composite modehng, implementing th e presc riptions 
of lSikora et alj (ll994l):'Dermer & Schlickeiser (2002'). can be 
found, for example, in Tramacere & Tosti (2003). The low 
state is described with the following "SSC+ERC" parameters 
B = 0.5 G, blob size i? = 7.9 x 10^^ cm, I? = 20, a log- 
parabola electron injection function between 7mi„ = 100 and 

Imax = 3 X 104, Ld,,k = 1-1 X TbLR = 0.1, T^sk = 

1.5 X 10^ °K and distance from the disk of 10^^ cm. The out- 
burst state is described with the following "SSC+ERC" pa- 
rameters B ^ 0.5 G, blob size i? 6.3 x 10^*^ cm, 2? 24, a 
log-parabola electron injection function between jmin = 100 
and jrnax = 3 X 10", Ldtsk = 1.1 x lO^^, tblr = 0.1, 
Tdisk = 1.5 X 10^ °K and distance from the disk of 7 x 10^'' 
cm. 

In Fig. \n\ the LAT averaged spectrum (E>200 MeV) of 
the entire Aug. -Dec. 2008 period considered (characterized 
by high statistics and consistent with an intrinsically curved 
shape described by a log-parabola model) is reported as well 
(the strip in the same orange color use d for non-simultaneous 
data). Log-parabola curvature (Section l24] i at high energy can 
be produced by several plausible models. For example by ra- 
diative particle cooling and stochastic acceleration processes 
driven by magnetic turbulence (rather th an systematic particle 
acceleration) acting near the shock front dStawarz & PetrosianI 
l2008h . If the energy where losses balance the acceleration 
rate if the acceleration time decreases more slowly than the 
loss time. The probability of energetic gain is lower when 
particle energy increases, because particles are confined by a 
magnetic field with a confinement efficiency decreasing for 
an increasing gyration radius. The integral energy distribu- 
tion of the accelerated particles results in a log-parabolic law: 
N{j) = A^o(7/7o) -P-i+''^ °g(T/^°), w here r is a curvature 
term dLandau et al. 1986; Fossati et al. 2000'; 'M assaro et al.l 
I2004t iPerlman et al. 2005; Tramacere et al. 2007h . Non- 
power law radiation spectra could also originate in a non- 
linear regime, such as in the presence of shock modifica- 
tions, precursors, a nd when diffusion coefficients vary with 
particle momentum (lAmato et al.ll2008]) . Another explanation 



takes into account episodic particle acceleration, that is appli- 
cable to high-energy flares with intrinsic spe ctral curvature 
( IPerlman et all 120051: IPerknan & Wilso"nll2005i) . The filling 
factor of the regions within which particles are accelerated 
is a function of both position and energy. If the light-crossing 
time of the emission region or the integration time of obser- 
vations is greater than the characteristic particle acceleration 
time, we effectively observe an electron distribution which is 
the product of a power law multiplied by a logarithmic term 
producing a spectral curvature. 

Another contribution to such deviations might also be the 
averaging over long periods of time, hence combining vari- 
ability effects, resulting in a cumulation of photons from 
different activity "flavors" (different source brightness and 
spectral hardness) being included in the evaluation of the 
spectrum. In the case of PKS 1502+106, the same spectral 
curved shape was also preferred for the much shorter, post- 
flare, interval, while the whole period includes many photons 
and covers a wider energy range (photons observed from the 
source have energies up to 15.8 GeV). PKS 1502+106, and 3C 
454.3 dAbdo et al.l2009ah for the first time show this departure 
from a simple power law, while EGRET observed only simple 
power laws, probably due to its lower high-energy observa- 
tional limit. The maximum peak observed in the LP models 
of PKS 1502+106 (in the i^F^y representation) are around en- 
ergies of about 390 MeV for the whole period and 480 MeV 
for the post-flare period. These peaks are consistent with a 
sub-GeV FSRQ blazar class. 

From Sect. 15.11 we remember that the inferred, apparent 
and isotropic, monochromatic luminosity at Eq — 100 MeV 
during the outburst phase (DoY 2008: 218.95 - 224.0) is 
Le>ioo — 1-1 X 10"*^ erg s^^, and the bolometric lumi- 
nosity is expected to be even higher than this value, mak- 
ing PKS 1502+106 one of the most powerful and luminous 
high-energy blazars observed during the first year the Fermi 
LAT all-sky survey. On the other hand The highest energy of 
photons detected from the source during the first four months 
of survey (15.8 GeV, assuming a strict PSF size criterium), 
has marginal consequences for extragalactic background light 
(EBL) predictions. The optical depth for 77 e+e^ pair 
production of 16 GeV photons propagating through the EBL 
from a redshift z — 1.839 source to Ear th approaches unity 
for rather high-density EBL models (e.g, IStecker et alj|2006l 
T^-y{z = 1.839, E =16 GcV) - 1.0 - 1.3), while most low 
density EBL models predict rather small interaction probabil- 
ities at such energies. An in-depth exploration of this finding 
will be presented elsewhere ( Abdo et al. 2009g). 

6. FINAL REMARKS 

This was the first time that PKS 1502+106, a distant radio- 
and X-ray-selected FSRQ, has been announced to have ob- 
servable high-energy gamma-ray emission above 100 MeV. 
Fermi LAT as basically a sort of all-space, -time, and -energy 
monitor allowed an ex celle nt spatial localization of this new 
7-ray source (Section 12.21 and FiglT] with a firm identifica- 
tion thanks to the optical-7-ray match of the outburst and, re- 
markably, also of a second flare). I t all owed detailed analy- 
sis of the energy spectrum (Section [2.41 and Fig. [TTl demon- 
strating the possibility for an intrinsic spectral curvature in 7- 
rays), and was able to pr ovide regular, daily monitoring flux 
light curve (Section [2. 3l l2]). This has made possible the dis- 
covery of consistent 7-ray brightness and activity ignited by 
the big outburst over the following 4 months (pointing out a 
variability behavior), the discovery of further minor 
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and rapid flares, and the disclosure of the outburst's tempo- 
ral shape. This type of "PSD-SED" monitoring performed by 
the LAT yielded advantage in developing the first unplanned 
Fermi multifrequency campaign, with a strategic 16 days of 
simultaneous Fermi-Swift monitoring. 

PKS 1502+106 is a powerful gamma-ray (~ lO"^ erg s'^ 
at i? > 100 MeV) FSRQ that showed, especially during the 
fast-rising outburst, a dominant MeV-GeV bolometric emis- 
sion similar to other FSRQs of the EGRET era. Dissipation 
probably occurred within the BLR, and, assuming for a black 
hole mass of ~ 10® Mq, the 7-ray emission was likely domi- 
nated by the ERC process. The SSC, in-jet, emission appears 
to dominate the observed SED from radio to X-rays bands. 
PKS 1 502+ 1 06 might be considered an example of a sub-Ge V 
peaked blazar, placed at the border of the BLR dissipated and 
the dusty torus/ambient-radiation dissipated FSRQs classes. 
The level of correlation found among the 7-ray, X-ray and 
optical-UV outburst and post-flare relaxing phases, support 
this idea. Opacity effects at cm and longer radio wavelengths, 
possible links between field ordering, jet-axis alignment, su- 
perluminal radio knots, and MeV-GeV outburst, are also de- 
picted by our results. 

In conclusion, the Fermi LAT performance in blazar sci- 
ence (as a stand alone observatory, or leading instrument for 
multifrequency campaigns), and the synergy between Fermi 
and Swift in particular, is evidenced by this work on PKS 
1502+106. By itself, this blazar is emerging as a major, lu- 
minous, energetic and 7-ray variable source, with promising 
diagnostic and discovery potential in emission modeling, in 
spectral and temporal variability studies, and in understand- 
ing the radio-gamma-ray connection. 
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